clear
delta=0.1;
n=0.02;
A=1;
k=[0.1;0.2];
sigma=0.1:0.01:0.9;
theta=0.36;

for j=1:length(sigma)
i=2;
while k(i)-k(i-1)>=0.01    
        k(i+1,1)=((1-delta)*k(i))/(1+n);
    i=i+1;
end
y1(j,1)=k(end).^theta;
end

% y=k.^theta;
% plot(k(1:299),k(2:300))
plot(sigma,y1)
